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jrt ' Abstract 

C/3 ' Based on the original idea of the density matrix renormalization group (DMRG) |l|, i.e. 

to include the missing boundary conditions between adjacent blocks of the blocked quantum 
Cy ■ system, we present a rigorous and nonperturbative mathematical formulation for the real- 

S ' space renormalization group (RG) idea invented by L.P. Kadanoff 12] and further developed 

by K.G. Wilson |H]. This is achieved by using additional Hilbert spaces called auxiliary spaces 
in the construction of each single isolated block, which is then named a superblock according 
to the original nomenclature |l| . On this superblock we define two maps called embedding and 
truncation for successively integrating out the small scale structure. Our method overcomes 
the known difficulties of the numerical DMRG, i.e. limitation to zero temperature and one 
space dimension. 



PACS: 75.10.Jm 



1 Introduction 

^^ ' Soon after K.G. Wilson's dramatic success in applying a momentum space formulation of the 

renormalization group (RG) method to the Theory of Critical Phenomena and the Kondo 
Problem B there was a considerable amount of efforts in applying the same type of approach as 
a real-space formulation to a variety of quantum physical problems. Since the momentum space 
'"^ ' formulation, apart from a few exceptions p, 0], relies in most cases on a perturbative expansion, 

C I real-space methods offer non perturbative approaches and are therefore extremely important in 

O ' applying RG ideas to strongly and complex correlated systems. It then turned out that for a 

variety of such physical models the real-space RG techniques give considerable bad results and the 
reason was unknown for nearly fifteen years. During that time some new real-space RG methods 
were discovered and some of them worked out very well whereas other methods failed without 

^ ■ giving any insight to their failure. We like to refer the interested reader to the book of T.W. 

5^ \ Burkhardt and J.M.J, van Leeuwen g] for a summary of work on this topic. 

Apart from these developments S.R. White and R.M. Noack published a series of papers containing 
a new idea for improving real-space RG techniques pi a. Based on the understanding of the 
importance of boundary conditions for isolated blocks in real-space RG methods for quantum 
physical systems a numerical approach was invented to take sufficiently many boundary conditions 
into account during the RG procedure. Apart from the impressive accuracy of the numerical results 
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this new approach displays also the typical universal character of a RG formulation, in that it is 

applicable with some particular changes for a variety of problems and was named the Density 

Matrix RG (DMRG) |l|, §. 

The dramatic success of the DMRG method has changed the picture of real-space RG techniques 

completely and has been applied until now in very different fields of scientific research |0, ^, ^■ 

The method itself is a rather complicated algorithm and a detailed description together with some 

examples is given by S.R. White [y. 

Despite of all the excitement concerning DMRG, the method has some important restrictions which 

are given by the method itself and therefore cannot be removed by applying simple changes to the 

DMRG algorithm. Here we mention the three main restrictions briefly: 

1. The chief limitation of DMRG is dimensionality Although higher dimensional variations are 
not forbidden in general, it becomes a complicated task. Recent applications of DMRG to 
finite width strips in two dimensions show a declining accuracy with the width. Therefore a 
successful approach for two dimensions in general or even higher dimension has never been 
worked out. 

2. DMRG is by definition an algorithm and therefore it is a purely numerical RG approach. 
Although this needs not to be a disadvantage we like to have an analytical formulation of 
the DMRG method. In such a reformulation the numerical DMRG scheme will occur as one 
possible realization of a more general description. We therefore expect a deeper insight to 
successful working RG approaches. 

3. DMRG is restricted to zero temperature and is usually applied for calculating ground state 
properties like the ground state magnetization or even the ground state itself. Finite temper- 
ature results were obtained only in the low lying spectrum but with very limited accuracy. 
In comparison to other real-space RG methods DMRG is different because it is designed 
to calculate ground state quantities. Recently, based on the idea of Xiang et al, a thermo- 
dynamic method was applied successfully, which combines White's DMRG ideaQ with the 
transfer-matrix technique[0 and which is now called TMRG. Although TMRG is also purely 
numerical since it shares the basic idea with DMRG it is an even more complicated algorithm 
||lC)| . Due to the close relationship to DMRG, the aim of TMRG is to give numerical accurate 
results for physical quantities and does not predict a RG flow-behaviour. In contrast our 
method is suited to calculate the flow-behaviour of the system, even analytically, although 
the main advantage by comparison with TMRG is the simple structure of our RG scheme. 
This makes it an easy task to apply it to a great variety of physical models. 

The rest of this article is organized as follows: In the next section we review the key idea of 
DMRG shortly. We begin by introducing the standard concepts of the real-space RG method 
in the language of spin chains in the way it was originally proposed. In section III we present a 
rigorous formulation of a real-space RG transformation. Each single block within the blocked chain 
is enlarged by an additional space, the auxiliary space. A single block together with its auxiliary 
space is called a superblock for which a real-space RG transformation is defined by integrating out 
the small spatial structure. Constructing a global RG transformation for the complete quantum 
system from concatenation of the local superblock RG transformations leads to the definition of 
exact and perfect RG transformations. In section V we give some final remarks including the 
relation to previous approaches in this direction. Applications in terms of this new formulation 
are shifted completely to a second paper. 



2 The idea of DMRG 

The very standard real-space RG approach is best explained for a spin Hamiltonian H on a one 
dimensional lattice as visualized in figure |l|. The dots represent the individual spins which are 
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Figure 1: A one dimensional spin chain. 

grouped together by breaking up the chain into blocks as visualized in figure |2| for a particular 
block composition of two sites. We like to establish a notation in which small letters refer to the 
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Figure 2: A one dimensional chain divided into blocks where each block is composed of two single 
sites. 



single site spins and capital letters denote the blocks. The block Hamiltonian for the block with 
the index / is then denoted as Hi . The idea of real-space RG is then to replace each block of the 
single spins by one effective block-spin, which leads to a renormalized block-spin Hamiltonian Hii . 
The calculation of the block-spins from the blocks composed of single site spins is carried out by 
a RG transformation R, which can be defined in various ways g, for example by projecting the 
block on the low lying spectrum [n] . In summary a RG approach is designed to split of the whole 
system into subsystems called blocks for which it is possible to reduce the degrees of freedom. 
Iterating this procedure leads to a RG flow in the parameter space of the model and the hope is 
to find a fixed point of this flow behaviour. Such a fixed point Hamiltonian is helpful to determine 
the universal behaviour of the physical model. 

As explained in the introduction, the boundary conditions of a block within the quantum system 
are essential for the calculation of a RG step, which is defined as one application of the RGT. In 
fact the different boundary conditions represent the correlations in the quantum chain between 
adjacent blocks to which we refer in the following as system blocks. To provide the opportunity to 
choose those boundary conditions, which result in the most accurate representation of an isolated 
block, the fundamental idea of DMRG is to embed the system block into a "bigger" block, called 
superblock. This nomenclature as well as the term system block is due to the original work of S.R. 
White ||l| . In some sense this simulates the environment represented by the surrounding spin sites 
and effectively smoothes out the sharp effects of the boundary conditions, as depicted in figure g. 
To construct a working approach out of this overall picture we are immediately faced with a 
twofold basic problem: How can we describe the embedding of the system block within the su- 
perblock and how can one include the boundary conditions during a RG step. In the framework of 
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Figure 3: An isolated system block embedded into a superblock visualized by dashed lines. 



DMRG these problems are overcome by focusing on one particular state, the target state \ijj), which 
is the ground state of the superblock Hamiltonian obtained by diagonalization. By using a complete 
set of eigenstates of the system block {|V'm*'*°™) i m — 1, . . . , ^system} and a complete set of eigen- 
states for the environment represented by the superblock { |^onvironmont^ ^ n—1,..., Zcnvironmont 
we decompose the target state {tp) according to 
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(1) 



We are interested in those states, which lead to an optimal representation of the target state \ip > 
in a "truncated" basis. Of course in this way we loose the exactness of relation (^ and we therefore 
denote the new result as an optimal approximation expressed as 
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where the optimal states {jV'p''*) i P—^i--- j ^opt < System} are defined in terms of the original 
system block states by 



IV'D = E - 



^opt 



Wl 



system 



) 



(3) 



with some coefficients aj^* . The coefficients 7p^„ in (0) can be determined by examining the 
reduced density matrix of the system block within the superblock [y . 

The twofold problem introduced above is therefore solved as follows: First the embedding of the 
system block within the superblock is achieved by reconstructing the target state of the superblock 
in a basis, in which the basis vectors are given as a tensor product composition of states of the 
system block and the chosen environment. In this way the system block is described within 
the bigger superblock. Since we have not truncated the set of the states which belong to the 
environment, the RG step for the system block is performed by taking all possible boundary 
conditions within the selected environment into account. 



From the previous discussion it becomes obvious that the coefficients 7p_„ can only be determined 
numericaUy within real applications of this technique. To develop a complete analytic approach, 
a method for using a target state will be impractical and we can only use the overall picture 
represented in figure & 

3 A rigorous real-space RG transformation 

We start this section by giving a very general but well known definition of a RG transformation 
(RGT). A ROT i? is a map defined on a set of physical variables {cr;} and a further set of parameters 
K = (ifi,X2,...) 

R: ({aJ,K)— .({a.„J,K') , (4) 

where {1} and {m} are not necessarily equal indexing sets and {nm} denotes the new set of blocked 
variables belonging to the larger scale. The quantitative prescription for the map (Q) is then given 
in physical terms by including physical constraints as for example the conservation of symmetries, 
the maintenance of the structure of the Lagrangian or the Hamiltonian, or the preservation of 
physical quantities, like for example the free energy of the system. Since in most cases it is a 
difficult task to define a transformation which combines all needed constraints this has led to an 
enormous variety of approximate RG transformations developed in the last decades |q] . 
The most common realization of the quantitative prescription is to apply the RG transformation 
R to the Lagrangian or the Hamiltonian as a functional which then acts on the variables and 
parameters given in (0). In the special example of a one dimensional quantum spin chain the 
new variables are the block spins and the new coupling constant belongs to the renormalized set of 
parameters K'. For our case we generalize this RGT to an arbitrary suitable functional dependence 
O 

R[0{{cJl},K)]=0{{^i^},K.') . (5) 

By further mathematical analysis of a particular RGT R defined by (ph this hopefully yields to a 
dependence of the renormalized parameters K' on the old parameters K which is called the flow 
behaviour of the RGT. We like to emphasize that once the functional dependence 0{{ai},'K) is 
known, we immediately know the functional dependence Cd/XmjjK') which plays an important 
role in our construction. 

We now make the Ansatz that in principle each RGT R can be written as a composition of two 
maps, called embedding and truncation ||l2]. This terminology originates from a RG technique for 



Hamiltonian systems [11 , which was then further developed and used for calculations of the flow 
behaviour |12| . Rephrasing equation (^ and focusing only on the renormalization of the set of 
parameters for determining the flow behaviour we get 

G+ o 0{K) o G = O (K') (6) 

where we denote G"*" as the truncation map and G as the embedding map. 

As a quite intuitive example for the abstract definition of the operators G'^ and G, in the special 
case of a functional dependence given by the Hamiltonian, we can construct G"*" as a projection 
map from the space of all eigenvectors of the Hamiltonian to a space containing a reduced number 
of eigenvectors. A projection map from this truncated space back to the space containing all 
eigenvectors is a natural way of defining G. Although such an example illustrates a possible 
application of the abstract formulation given by (g) it raises the question of which eigenstates are 
necessary to keep for constructing the truncated space. In the case of zero temperature we can 
argue that only those eigenvectors should be kept, which correspond to the low energy eigenvalues 
[g . As pointed out previously the idea of this paper should be to invent a real-space RG formulation 



which overcomes these hmitations by a more abstract formulation. 

Let us now assume that the functional dependence O is given by some operator, not necessarily the 

Hamiltonian, on the original Hilbertspace Ti, so that equation (o) can be written as the commuting 

diagram 

n' ^ ) n 

0(K')| |o(K) (7) 

n' < n 

G+ 

where TC' refers to the effective Hilbertspace for the functional dependence of the truncated set of 
parameters. 

We introduce the blocking concept discussed in the previous section as a tensor product decompo- 
sition of the Hilbert space 

n = (^Hi, (8) 

led 

where 3 denotes some indexing set for the blocks. We are looking for an embedding and truncation 
map which respects the block decomposition by factorization 

Gw = (S>Gn', and G+ =(S)Gtir- (9) 

led ie3 

Using this mathematical formulation of the blocking scheme we like to write the RG transformation 
for a block in an analogous way 

Oh'^ (K') = G+^ o On, (K) o G«, , (10) 

due to (p) . But equation (HG) is not an independent relation since we have to relate it to the global 
relation (ph. To decompose (P) into the blocked pieces (10) we have to assume that the operator 



O-H can be decomposed into commuting block operators O-j-c, which is not the case in general in 
quantum physics. Therefore the problem encountered so far is to find suitable functions Ou (K) 
which respect the block decomposition of the Hilbert space within the RGT. 

To find a solution for this problem our Ansatz is to enlarge the Hilbert space Ti by an additional 
(auxiliary) Hilbert space Tiaux due to the composition rule 

Htotal ^n® Haux . (11) 

We like to think of the space Titotai as some kind of 'super space' and the global operator 
Cw^Waux (I^) is then embedded into the total space Titotai- The key idea is to recover a block 
decomposition for O-H^Waux C^) i^^o blocked pieces of the form ^'^^^(•h^^^) (K) which we identify 
as superblocks according to section 0. The next step in our approach following the basic principles 
of DMRG is to outline a general construction for Onign^^^ (K) with a commuting block decompo- 
sition. This can be performed explicitly by starting with standard real-space RG concepts. 
In the formulation of standard block RG we consider a decomposition of On^n^^^ (K) into discon- 
nected block functions given by 

0«,«(«_)^(K) = o;,^;;^-_)^(K) with /e a. (12) 

where we have neglected the non commutativity or correlations between the blocks completely. A 
straight forward way to improve the standard RG method is to include somehow the correlations 
between adjacent system blocks. As visualized in figure H we can refer to these correlations as 
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Figure 4: Successive blocks in a one dimensional chain in the commuting case a) and the non 
commuting case b). In the non commuting case the system blocks are connected by correlation 
blocks drawn by a dashed line. 



blocks, which we denote therefore as correlation blocks. Using these correlation blocks enables us 
to represent the non commutativities between the system blocks in a compact way and we denote 
these correlation blocks, using the overall notation given in the appendix, as 



O: 



correlation 
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(K) 



(13) 



with U{i.i-i^... } «) (Haux){i,i_i,... } C Ti/ «) {Hi,ux)i ® H/_l (Haux)/_i 

The subspace 

W{i,i_i,... } «) (Waux){j^i_i^... } = Hi ® (Waux)j «> Hi^l «> (Haux)j_i ® • 



(14) 



denotes the tensor product composition of all the block Hilbcrt spaces used for the construction of 
the correlation block. 



4 Decomposition rules 

We are now dealing with the problem how to include these correlation blocks into the RG trans- 
formation. One can find approaches in the past where this is performed perturbatively [|l2[ and 
therefore unsuitable in our case. To find some insight into this problem let us start with the 
composition 
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(15) 



which is exact and as always {i,i — 1, . . .} denotes the subset of ah needed product subspaces 
for constructing the correlation blocks. We like to stress that the decomposition dl3) in sums of 
system blocks and correlation blocks is not unique. Later we will examine another decomposition 
to which in contrast we will refer to as the product decomposition. 
Let us apply the RG transformation (0) on the sum decomposition (p^) 
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(16) 



and all quantities are used in the context of additional auxiliary spaces. Let us first consider the 
system block summand in (|lq), which can be rewritten as 
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(17) 






This is exactly the local RGT for the system blocks (10) if we neglect the last product term on 
the right hand side of (|lj). We will refer to this factor as a correction term that vanishes if we 
demand 
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(18) 



Inserting this constraint into (O) , carrying out the same calculation for the correlation blocks and 
finally using relation (O) we get the renormalized version of equation (ITq) given by 



ou'^w^^^ (K') = Y. o^-^^;f^,_^^ (K') + Y. o«"tr°"^ 
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which leads to the renormahzed set of parameters K'. In (|T^) we used the reasonable definition 
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Relation ( [iq ) introduces an additional constraint for the RGT and therefore restricts the variety 
of possible transformations. 



In the case of a product decomposition of the operator 0(K) we can write 
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In analogy to the case of the sum decomposition (|T^) we can apply the RG transformation (g) to 
( pl| ) which leads to the expression 
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Since this is already the final step in the calculation for this special case of a decomposition we are 
not able to write the result as a composition of the renormalized system block part and correlation 
block part as we did in ( ]l9| ) for the sum decomposition. By the considerations so far the product 
decomposition therefore seems to be not as useful as the sum decomposition for later applications. 
This is not the case as we will show in the following. 

For the auxiliary space we distinguish between two different cases, an active role and a passive 
role. Here active means that the auxiliary space is directly involved into the RGT, i.e. G and G^ 
act nontrivial on this additional space. The commutative diagram describing the general active 
situation is given in (123 
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(23) 



Relation (p3) reduces to a rewriting of (R) , if the transformation maps G and G+ each operate as 
the identity on the auxiliary space and the functional dependence O (K) acts non trivial only on 
7i. This gives us an example of the particular case of a passive role of the auxiliary space as it is 
depicted in (24). 



/®i^ 



H'®Kux 






ux 


n ® Haux 












0>l'»«Ux(^') 








OwgHau 


x(K) 


'H <S> Haux 








n ® Haux 




'■ 










G+ 


®l«ai 


X 







(24) 



In the case of (ESh we can think of the auxiliary space as some kind of medium not changed during 
a RG step. The active and the passive choice of the auxiliary space yield two different realizations 
of our RG, which we will call the 'general (real-space) RG' (GRG) and refer to the corresponding 
RG transformation as GRGT. 



5 The construction of the local GRG transformation 

So far we have discussed different types of quantum decompositions and types of auxiliary spaces. 
We now turn to the question how to construct the embedding map Gu'^w ^^nd the truncation 
map G^^^ ■ In (|^) we used the functional dependence O to introduce physical constraints within 
the RG transformation. To determine G-h'^w and Gj^^^ we introduce another constraint. 
In addition to keeping the structure of the operator O we relate O to a physical quantity Z{0) 
which acts as a physical invariantn. Equating the original physical quantity Z{0) calculated from 
the original quantum lattice and the effective physical quantity Z{0') for the reduced lattice we 
obtain G^^^ and G-h'^h'^^^ from 

Z [0«®w_ (K)] = Z [G+^„_ o 0„«H_ (K) o G«««_] = Z [On'^m^^ (K')] . (25) 



We refer to equation (ES^ as the invariance relation for the RGT. Finally we have to decompose 

'^ti'sn^^^ and Gu'®n'^^^ according to (|). 

We are now able to give the precise definition of the local RGT in the form 



'^«/<S(-«aux)i 



n'j ® (Kux)/ °" ' ni ® in.^^)j 



°«i8>(«Ux).(^')| |o„^«(„._)^(K) (26) 



K ® (Kux)/ -^ ni ® (Waux), 

where we refer to Gt, „ ,^ -. and G^' ^ mi \ as the generators of the transformation. By the 
explanations of section ^ 

0«.«(«.„x). (K) === OCT«_)^ (K) • n ^«T:,t";^(«..x),......>(K) (27) 

CI 

or 0«,«(K„x). (K) = 0;,^;7^_)^ (K) + J2 ^«T:,t";«(«au.),..,..->(^) (28) 

CI 

and analogously for C7^^®(w^„j, (K'). 

6 Perfect and exact local RG transformations 

In this section we study the relationship between (|2^) and the global RGT 

n' ® K^, — ^^^ n ® Tiaux 

'^«'»«Ux(^')| |o«®(Ha„x)(K) (29) 

^'•^Kux * H®Haux 

Diagram (p9) represents an exact relation which implies all the necessary constraints for the RG 
procedure as can be verified from equation (|2^). We therefore choose relation (g^) as the basic 

t A possible example for such a quantity can be the partition function or the free energy of the physical system. 
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relation in defining local RGTs. 

Decomposing the global RGT (E9h into local RGTs given by (Eff) demands for a decomposition of 
O into commuting blocks. From previous considerations we conclude that this is impossible for 
quantum chains due to the correlation blocks occurring in a decomposition of a quantum physical 
system. Therefore the idea is to use the auxiliary space to decompose the chain into commuting 
blocks by storing the information about the correlations of adjacent system blocks within the 
auxiliary space. By the decompositions discussed so far we then decompose a chain into system 
blocks and try to find an auxiliary space (Tiaux)/ for each system block which takes over the role 
of the correlation blocks within the RGT as visualized in figure |^. This statement can be made 
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Figure 5: A rigorous blocked chain in the non commuting case with a) decomposition into system 
and correlation blocks, the latter visualized by dashed boxes and b) decomposition into system blocks, 
each equipped with an auxiliary space suitable to take over the role of the correlation blocks during 
the RGT. 

more precise. 

Definition 6.1 A local RGT is said to be perfect if there exists a local operator 



'^'H'M'H'au.) I (^') = 






° O'l^'ZTn...)^'^) ° 



G-H'j ® (n'^,J, 



together with a global functional dependence O-u'^n' (K') defined by the decomposition 

Ou'^w^^^ (K') := Y. On'MUL.), (K') or dw^w^^^ (K') := [] On'Mn'.^J, (^0 
and no further local relation governing the renormalization of the correlation block part occurs. 
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The main advantage of a perfect RGT is a rigorous mathematical description for a local RGT. 
Although the structure of the local Operator Oui®(h.^^^) is conserved, a perfect RGT does not 
make use of the invariance relation ([25|). 

Definition 6.2 A local RGT is said to be exact if it is perfect and 

Z [Ow®H„„ (K)] = Z [Ow'87Y'_ (K')] = Z [Ow'®w'_ (K' 

If a RGT is exact it includes all needed constraints and therefore we can compare the RGT to the 
classical situation where non commutativity effects are absent. 

At this point we like to give some important remarks on perfect and exact RGTs. Although in 
both cases a rigorous mathematical formalism is used, a physical approximation usually enters the 
problem by choosing an appropriate auxiliary space. Only for a certain class of models we will 
be able to find auxiliary spaces with a structure that allows for describing the non commutativity 
effects without any approximation. 

We stress again that in the exact as well as in the perfect RGT O-h and O-h' are known so that 
we can determine G and G^ in both cases according to the explanations in section |5|. 
If the auxiliary space is active it may happen that it vanishes by truncation during the RG proce- 
dure. In such a case no auxiliary space is available after the local transformation has been worked 
out and the previously provided information concerning the correlations between adjacent system 
blocks is lost. Therefore the RGT is at most perfect. 

In the case of an auxiliary space which (only) allows for an approximate description of the cor- 
relations between adjacent system blocks we would like to have some insight into the accuracy of 
the approximation. Here we remember the numerical DMRG procedure in which convergence of 
numerical values of ground state quantities by enlarging the superblock is used as an estimate for 
the accuracy of the method. 

It is apparent that only in the case of an exact RGT we are able to calculate global quantities like 
the total ground state energy shift. Since we are mainly interested in an overall effective coupling 
determining the RG flow we are looking for exact RGTs. 

7 Conclusions 

We have invented a non perturbative quantum RG method based on the idea of an additional aux- 
iliary space. The work was motivated by the success of the DMRG concerning numerical results 
and the open question of an underlying general mathematical framework. 

The main objects introduced in this article are the auxiliary space Tiaux and the two maps 
G^ « (n ) ^^d ^H' «) {W^^J which generate the RGT. By using these quantities we were able 
to give the definition of an exact local RGT which is the final result of this work. An exact local 
RGT involves all the information provided by the physical system. 

In future work we will proceed by applying our abstract formalism presented here to quantum 
spin chains like the Heisenberg models and compare our results in the context of related work on 
these models |0|. This leads us to concrete and different examples of possible auxiliary spaces. As 
expected, the correct choice of the auxiliary space will be the main ingredient in the construction 
of the RGT, whereas the definition of the maps G^ „ ,^ > and G-^' ® (w^^^) turns out to be 
rather straight forward. We also hope for further applications of the method introduced here. 
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Appendix 

Throughout this work blocks are denoted by capital indexing letters, corresponding to the block 
sites. The indexing set for the blocks is denoted as 3. Neighbouring blocks are denoted by a se- 
quence /, / — 1, / — 2, • • • G J whereas arbitrary blocks are indexed by different letters I, J, ■ ■ ■ € 3. 
A block Hilbert space Tii contains at minimum two single site Hilbert spaces Tii and Tii-i- Single 
site Hilbert spaces are denoted by letters i, j, fc, . . . . To point out that a single site space TLi is 
contained in a block space TLi we write Tii C Tii or even simpler i G / if it is clear that / refers 
to the block Hilbert space. We also use the abbreviation {i, i — 1, . . . } C {/, / — 1, . . . } instead of 
writing 'Hu,i-i,,,, \ C "Hsn ® Tiij^iy, .... By this notation it becomes not clear which single site 
space is contained in a certain block Hilbert space. If this is important it must be pointed out 
explicitly. 

Expressions which are written in the form expression arc either defined and used in this work or 
have special physical meaning. 
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